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ABSTRACT 

Rex1/Zfp42 is a Vyf-related zinc-finger protein 
whose expression is frequently used to identify 
pluripotent stem cells. We show that depletion of 
Rexl levels notably affected self-renewal of mouse 
embryonic stem (ES) cells in clonal assays, in the 
absence of evident differences in expression of 
marker genes for pluripotency or differentiation. By 
contrast, marked differences in expression of 
several endogenous retroviral elements (ERVs) 
were evident upon Rex1 depletion. We demonstrate 
association of REX1 to specific elements in 
chromatin-immunoprecipitation assays, most 
strongly to muERV-L and to a lower extent to IAP 
and musD elements. Rex1 regulates muERV-L ex- 
pression in vivo, as we show altered levels upon 
transient gain-and-loss of Rex1 function in pre- 
implantation embryos. We also find REX1 can asso- 
ciate with the lysine-demethylase LSD1/KDM1A, 
suggesting they act in concert. Similar to REX1 
binding to retrotransposable elements (REs) in ES 
cells, we also detected binding of the REX1 related 
proteins YY1 and YY2 to REs, although the binding 
preferences of the two proteins were slightly differ- 
ent. Altogether, we show that Rex1 regulates ERV 
expression in mouse ES cells and during pre- 
implantation development and suggest that Rex1 
and its relatives have evolved as regulators of 
endogenous retroviral transcription. 

INTRODUCTION 

After undergoing a first differentiation step, the pre- 
implantation blastocyst is divided in Inner Cell Mass 
(ICM) that gives rise to the embryo proper and 



trophectoderm, an external epithelium that contributes 
to the placenta. Self-renewing stem cells that can be 
derived from each of these lineages (1) are referred to as 
embryonic (ES) and trophectoderm (TS) stem cells, re- 
spectively. ES cells can be maintained in culture for an 
apparent unlimited number of cell divisions (self-renewal) 
and maintain the defining property of pluripotency or the 
ability to differentiate into cell lineages of all three 
primary layers of the embryo. Molecular mechanisms 
that maintain this pluripotent self-renewing state operate 
at different levels and include (but are not limited to) 
signalling by leukaemia inhibitory factor (LIF) and 
BMP4, inhibition of ERK signalling (2), co-operating 
networks of transcription factors and epigenetic mechan- 
isms (3,4). The Oct4, Sox2 and Nanog transcription 
factors (5-8) constitute a core transcriptional network to 
maintain pluripotency through mutual positive regulation 
(4) and collaborative regulation of target genes. A distinct 
module whose function is essential for the maintenance of 
pluripotency and self-renewal consists of Cnot3, Trim28, 
c-Myc and Zfx (9,10). Trim28 has recently also been 
shown to participate in the repression of endogenous 
retroviral elements (ERVs) in mouse ES cells (11). 

Rexl was first discovered as a result of its specific ex- 
pression in pluripotent F9 embryonal carcinoma (EC) 
cells (12). Rexl (for reduced expression- 1, also known as 
Zfp42) was subsequently shown to be expressed in other 
pluripotent cell types, especially undifferentiated ES cells 
(13), multipotent adult progenitor cells (14) and amniotic 
fluid cells (15), in the germ cells of the testis and in the 
ICM and early trophectoderm derivatives of the mouse 
embryo (13). Rexl expression has been positively linked 
to increased pluripotency in both mES cells (16-18) and 
human ES and iPS cells (19,20). In contrast, conflicting 
results have been reported regarding the functional role of 
Rex-1. Gene silencing by RNA interference results in loss 
of self-renewal in ES cells (21) and overexpression of Rexl 
negatively affects self-renewal (D. Guallar, M. Sanchez 
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and J. Schoorlemmer, unpublished data). However, Rexl 
does not have to be provided for efficient reprogramming 
of differentiated cells towards iPSs (16,17), Rex-1 is dis- 
pensable for maintenance of self-renewing pluripotent ES 
cells (22) and ES cell lines can efficiently be derived from 
itex7 -deficient blastocyst (23). 

Rexl encodes a protein containing four Cys-His type 
zinc fingers, which is localized in the nucleus in ES cells 
(23), and displays significant similarity to the YY1 tran- 
scription factor family in the DNA-binding zinc-finger 
domains (24). Rexl target genes have been surveyed by 
gene association and differential expression studies. 
Target genes identified in ES cells now encompass a 
circuit of active genes implicated in protein metabolism 
that coincides partially with Myc targets as opposed to 
Oct4/Nanog/Sox2 targets (25), binding to (and regulation 
of) Tsix regulatory elements (26), RinglBj Rnf2 regulated 
genes (27) and also imprinted genes during pre- 
implantation development (28). Interestingly, the 
absence of Rexl from an ES cell line has been linked to 
loss of pluripotency upon prolonged passage and to 
enhanced expression of retrotransposable elements 
(RE) (29). 

Transposable elements (TEs) are repetitive DNA se- 
quences that are ubiquitous and abundant components 
of most genomes including mammals and constitute 
>45% of the human and mouse genome (30,31). They 
can duplicate and reinsert within genomes, either autono- 
mously or with the assistance of proteins encoded by other 
(related) elements. As a result, TE profoundly impact 
genome function and evolution, as transposon-derived 
promoters also direct expression of alternative transcripts. 
De-regulated gene expression mediated by the activation 
of transposon promoters contributes to tumorigenesis and 
autoimmune disease (32,33). 

Most mammalian TEs are REs, which propagate 
through an RNA intermediate, and 8-10% of those are 
retrovirus-like long terminal repeat (LTR) elements 
referred to as ERVs. They constitute a range of similar 
but clearly distinguishable elements with varying copy 
numbers, autonomy and expression patterns that 
together occupy ~5.4% of the mouse genome (34). The 
superfamily of ERVs is composed of (but not limited to) 
muERV-L, IAP, musD, ORR1 and MT families 
(Supplementary Table SIV), with varying copy numbers 
ranging from 300 to 200 000 (31,35). Transcription of 
ERVs in different species is elevated in germ line cells, 
early embryo and placenta compared with adult or 
differentiated tissues (36). MaLR and muERV-L family 
ERV are frequently encountered in chimeric transcripts 
in the mature oocyte/zygote and two-cell embryos, re- 
spectively (37). MuERV-L also displays a transient 
increased expression during the two-cell stage (38) and 
normal development to four-cell embryos is compromised 
when muERV-L levels are attenuated (39). By contrast, 
IAP and musD/ETn elements are specifically expressed 
from the blastocyst stage onward (40). As expression of 
several families of elements is subject to tight regulation 
during embryonic development and each element was 
found to act as promoter for subsets of genes at defined 



stages (37), it has been proposed that MaLR, muERV-L 
and IAP synchronize stage-specific gene expression. 

ERV transcription is limited in most tissues due to 
methylation-directed silencing. Pre-implantation develop- 
ment is accompanied by general demethylation (41), 
calling for specific mechanisms to regulate expression of 
ERV elements that include histone modifications 
[reviewed in (42)] and RNA interference (43), protection 
from demethylation (44) or specific repressors including 
KAPl/Trim28 (11), ESet (45) and LSD1/KDM1A (46). 
Some of these mechanisms that constrain expression 
during pre-implantation development are also operative 
in ES cells (11,45). In contrast to 762/?/-deficient ES cells 
that lose pluripotency in a few generations, Lsdl -deficient 
ES cells can be stably propagated although muERV-L is 
severely de-repressed (46). 

Expression of Rexl in the germline (12,47), placenta 
(28) and during pre-implantation development (13) coin- 
cides with high prevalence of gene expression directed by 
retroviral elements. The presence of Rexl in tissues dis- 
playing high ERV expression prompted us to consider a 
potential direct role of Rexl in regulation of REs. 

MATERIALS AND METHODS 

Immunological reagents and western blot 

The oc-YY2 and the oc-REXl sera raised in rabbit have been 
described previously (27). The amino-terminal region of 
REX1 (amino acids 1-116) that we have used for immun- 
ization does not show homology at all to the corresponding 
region in YY1 (Supplementary Figure SI A) and does not 
cross-react with YY1 (Supplementary Figure SIB). The 
rabbit oc-REXl serum was further affinity-purified over 
REX1A-GST protein (ImmunoStep SA, Salamanca, 
Spain). Monoclonal oc-HA (clone HA-7) was obtained as 
an unpurified ascites fraction (Sigma H9658). 

For western blot analysis, cells were scraped with 
pre-cooled phosphate buffered saline (PBS), pelleted by 
centrifugation and re-suspended in lysis buffer [50 mM 
Tris (pH 8.1), 10 mM EDTA, 1% (w/v) sodium dodecyl 
sulphate (SDS)] containing protease inhibitors (PMSF 
and Complete; ROCHE). To solubilize REX1, samples 
were sonicated (Bioruptor®, Diagenode) until homogen- 
ous in ice water (0°C). Debris was removed by centrifuga- 
tion for 15min at 13 000rpm at 4°C. Laemmli sample 
buffer was added to lysates, and samples were separated 
by SDS-polyacrylamide gel electrophoresis and analysed 
by standard western blot as described in Supplementary 
Materials and Methods section. 

Co-immunoprecipitation between REX1 and LSD1 was 
assayed as described in Supplementary Materials and 
Methods section, using extracts from 293T cells 
co-transfected with either HA-tagged REX1 or with a 
mixture of plasmids expressing HA-tagged REX1 and a 
Flag-tagged LSD1A construct (48). 

Plasmid construction 

For the generation of short hairpin RNA (shRNA) 
vectors, the following fragment were fused into the 
BamHI-Bglll backbone derived from PGKHygro 
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(pHPCAGGS): a doxycyclin-inducible HI promoter (49), 
shRNA as Bglll-HIII fragment, and a promoterless eGFP 
derived from pCH-Octi (50). Constructs were verified by 
DNA sequencing. The shRNAs were designed as 
described (51), and sequences are as follows (shlRexl: 
5'-ACGGATACCTAGAGTGCATCA, sh2Rexl: 5-CA 
CGGAGAGCTCGAAACTAAA, shRNA Gfp: 5'-AAG 
CGCGATC AC ATGGTCCTG) . Plasmids used for trans- 
fections were purified on PureLink™ kits and columns 
(Invitrogen). 

For overexpression, a Rexl/Zfp42 complementary 
DNA (cDNA) (27) was fused to an IRES-eGFP derived 
from pIRES2-EGFP (Clontech) and inserted into the 
chicken ft Ac tin promoter-driven expression vector 
pCAGIP (52). Further details of all plasmids used 
are available upon request. Plasmids used for injections 
were linearized and recovered from GeneClean Turbo 
cartridges (MP Biomedicals) in Tris-EDTA (10/0.1 mM) 
as described (M. Climent et ai, submitted for 
publication). 



Cell culture, differentiation and transfection 

293T cells were maintained and transfected with mixtures 
of plasmids according to the manufacturer's protocol 
(Lipofectamine 2000; Invitrogen). 

Mouse TS cell line B7 (53) was maintained as described 
(54) on 1% gelatine-coated tissue culture dishes supple- 
mented with FGF4 (25 ng/ml)(Peprotech) and heparin 
(1 mg/ml)(Sigma). ES cell line E14T (55) was maintained 
on gelatine-coated tissue culture dishes in medium supple- 
mented with 10% foetal calf serum, LIF and 2i as 
described (56). To block HDACs or DNA methylation, 
ES cells were cultured for 24 h in the presence of 200 nM 
trichostatin A (TSA; diluted from a 1 mM stock solution) 
or 3 |iM 5-aza-cytidine (diluted from freshly thawed 
aliquot). RA (Sigma) was stored as a lOmM stock in 
dimethyl sulfoxide (DMSO) at -80°C, freshly diluted 
before use at 0.1 mM in culture medium, and added 
directly onto the cells to a final concentration of 1 uM. 
Cells were aggregated in the presence of DMSO as 
described (57). E14T ES cells were transfected by 
standard electroporation (5 x 10 6 cells; 20 ug DNA; 
200 V, 960 uF, ooQ). ES cells electroporated with 
plasmids expressing shRNAs were seeded at clonal 
densities in the presence of LIF and 2i to allow for 
colony formation as a measure of self-renewal (usually 
100 000 cells per six well). After 24 h, transfected cells 
were subjected to selection (H3274; Sigma at 160ug/ml) 
for the Hygro resistance carried on the same plasmid. 
After selection for 7 days, cells were fixed overnight at 
4°C in 80% ethanol, and stained for alkaline phosphatase 
(AP) using a Sigma kit (85L3R). Alternatively, cells were 
harvested and processed for expression analysis or chro- 
matin immunoprecipitation (ChIP). AP+ colonies were 
counted in fields of predetermined size on a Nikon 
Diaphot Inverted Tissue Culture Microscope. For each 
condition, 10 fields were counted, colony numbers are rep- 
resented as mean ± standard deviation (SD). 



RT-qPCR analysis 

Cells were washed with PBS, scraped and total RNA was 
extracted using TRIzol® reagent (Invitrogen). After diges- 
tion of genomic DNA (RQ1 RNAse-Free DNase; 
Promega), RNA was extracted with phenol/chloroform, 
precipitated with ethanol, re-suspended in water and 
quantified using Nanodrop (Thermo Scientific). cDNA 
was synthesized from 2 jig of RNA either with Oligo dT 
or random hexamer primers (Thermo Script® RT-PCR 
System; Invitrogen) and stored at — 20°C until used. 
cDNA was analysed by quantitative PCR (Platinum® 
SYBR® Green qPCR SuperMix-UDG; Invitrogen) on 
an ABI Prism 7000 Real-Time PCR system, reactions 
were performed in triplicate. Data were processed 
using the AAC t method (58), using Gapdh as a reference 
gene. Sequences of all primers used are included in 
Supplementary Materials and Methods section and 
Supplementary Table SI. An extensive justification of 
the localization and origin of primers used for 
RT-qPCR of TE is provided in Supplementary 
Materials and Methods section. 

siRNA used for microinjection 

Stealth siRNA duplexes against mRexl and control 
siRNAs (Ref 12935-112) were purchased, annealed ac- 
cording to the manufacturer's instructions (Invitrogen) 
and changed to lOmM Tris/EDTA 0.1 mM as described 
in the accompanying paper (M. Climent et ai, submitted 
for publication). 

Oocyte and embryo collection, culture, microinjection and 
RNA isolation 

All procedures were carried out under Project Licence 
PI29/08 approved by the in-house Ethic Committee for 
Animal Experiments from the University of Zaragoza 
(Spain). 

Oocyte and embryo collection, culture and microinjec- 
tion were performed according to standard procedures 
(59). A detailed description of microinjection procedures, 
embryo culture, Rexl overexpression and gene-expression 
analysis in embryos will be provided elsewhere (M. 
Climent et aL, submitted for publication). One-cell 
embryos were microinjected with siRNAs at 100 uM or 
linearized plasmid at 2ng/jil and 0,4mg/ml Dextran- 
coupled Texas Red (D1829; Invitrogen). In brief, 
embryos were injected on Day 1 , selected for development 
to at least the two-cell stage on Day 2 and to the five- to 
eight-cell stage on Day 3. Embryos injected with pCAG- 
itex7IRESeGFP were separated on Day 3 in three groups 
representing low (absent), intermediate and high eGFP 
expression. Only the groups identified as low and high 
expressors were used in the experiments. 

For the analysis of mouse pre-implantation embryos, an 
identical number of experimental and control embryos 
from the same experiment was used (typically 8-10 
embryos). Total RNA was reverse-transcribed and two 
oocyte or embryo equivalents were used as a template 
for each PCR. PCRs were carried out for 35 cycles of 
denaturation at 94° C for 30 s, followed by 30 s of 
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annealing, 1 min extension at 72°C; final extension 5min. 
Products were separated on agarose gels, visualized using 
ethidium bromide and photographed on a Gel Doc 
transiluminator (BioRad). Data were quantified using 
Quantity One software (BioRad), and expression levels 
were recalculated using H2afz as a reference gene. 
Primers used in this study are shown in Supplementary 
Materials and Methods section and Supplementary 
Table SI. 

ChIP and locus-specific PCR 

ChIP assays were carried out as described (27), using per 
5 x 10 6 cells the following sera for immunoprecipitation: 
Pre-immune serum (75 jig), rabbit anti-REXl IgG (75 jig), 
affinity purified rabbit anti-REXl (25 jug), mouse mono- 
clonal H-10 anti-YYl (Santa Cruz; 2.5 ug), rabbit 
anti-YY2 (25 ug). We compared semi-quantitative PCR 
amplification on amounts of chromatin obtained from 
the same number of cells after immuno-precipitation 
using either pre-immune serum or ocREXl serum. PCR 
products were visualized using ethidium bromide and 
photographed. Primers used are listed in Supplementary 
Materials and Methods section and Supplementary 
Table SIL 

Quantitative real-time PCR on ChlPs was performed in 
96-well plates using Platinum® SYBR® Green qPCR 
SuperMix-UDG (Invitrogen) on an ABI Prism 7000 
Real-Time PCR system and processed using the AC t 
method. Amplification was performed in 15ul reactions 
using the following parameters: 50°C for 2 min, 95°C for 
lOmin. followed by 40 cycles of (95°C 15s, 60°C 60s). 
Data presented are derived of triplicate reactions; the 
results shown are the aggregate of a minimum of three 
separate ChlPs performed on different days. Data were 
processed to calculate percentage association as described 
(www.SABiosciences.com) and are represented as the 
mean ( ± SD/SEM). 

Real-time PCR analysis 

Standard curves of all primers (used at 200 nM) were per- 
formed to check for efficient amplification (between 90 
and 110%). Melting curves were also performed to 
verify production of single DNA species with each 
primer pair, except for ERV sequences. Relative levels of 
expression in each assay were obtained through the 
standardized AAC t method. We used (i) Gapdh mRNA 
levels as a reporter in all experiments of Rexl knockdowns 
and (ii) the appropriate control cell line (DMSO-treated 
aggregates) as the reference sample in Figure 2B. Data 
presented are the mean of three to five independent experi- 
ments each performed on different days and using a dif- 
ferent preparation of cells. The relative amount of 
expression between samples was compared using the 
Tukey's method (95% confidence interval) following 
one-way ANOVA (Instat version 3.0, GraphPad 
Software, San Diego, CA, USA). 

Enrichment levels in ChIP assays are expressed as a 
percentage of immuno-precipitation relative to the input. 
Essentially, the AC t method was used to calculate a ChIP 
over input ratio, after subtraction of C t s derived from 



control immunoprecipitations using pre-immune serum. 
Values obtained were corrected for the appropriate 
dilution factor of each analysed fraction, and multiplied 
by 100 to get the percentage of immunoprecipitation 
(% IP). Alternatively, data are shown as fold enrichment 
of Rexl binding for a particular locus compared to control 
loci as indicated. 

RESULTS 

Depletion of Rexl expression causes a loss of 
self-renewal in mouse ES cells 

Although it has been reported that Rexl ( Zfp42 ) null mice 
appear normal, viable and fertile and normal ES cells can 
readily be derived from Rexl null mice or by homologous 
recombination in vitro (22,23), Rexl expression is 
associated with highly pluripotent populations of stem 
cells (18,20), suggesting an important role for Rexl in 
ES cells. We designed shRNAs to attenuate Rexl expres- 
sion, which were cloned in appropriate expression vectors. 
Resulting vectors were tested by co-transfection with 
HA-Rexl in 293T cells. Results show that expression of 
HA-REX1 (Figure 1A, lanes 1 or 4) was abolished by 
co-transfection of either of the two hairpin constructs 
(termed Shi Rexl and ShlRexl). Transfection efficiency 
was measured using a co-transfected eGFP plasmid and 
was roughly equal among samples (23-37%). 

We subsequently investigated Rexl function in ES cells 
by using stable episomal expression of these shRNAs. ES 
cells were electroporated and selected for expression of 
shRNAs, seeded at clonal densities and colony formation 
was counted as a measure of self-renewal. As expected, no 
colonies were observed after selection of cells 
mock-transfected in the absence of plasmids that carry 
the HygroR gene (data not shown). ES cells that express 
shRNAs targeting Rexl (Figure IB) clearly formed fewer 
colonies when compared with cells electroporated with 
control shRNAs, consistently giving rise to 44 and 67% 
of colonies in control assays for Shi Rexl and ShlRexl, 
respectively (Figure IB and C). All colonies remained un- 
differentiated as judged by morphology and AP staining 
(Figure IB, top panels). The loss of colony formation was 
not caused by off-target effects, as Shi Rexl did not nega- 
tively affect colony formation in Rexl deficient cells 
(Supplementary Figure S2). 

We did not observe complete elimination of Rexl mRNA 
levels, but a clear reduction to 33% or 45 % of wild- type 
levels for Shi Rexl or Sh2Rexl -expressing ES cells, respect- 
ively (Figure ID). More important, REX1 protein levels in 
E14T ES cells (Figure IE, lane None) were severely reduced 
in cultures overexpressing either hairpin RNA (Figure IE, 
lanes marked ShRexl), as opposed to cultures expressing 
control ShRNA (Lane Con). Together, these data show 
that reduced expression of Rexl was coincident with dimin- 
ished self-renewal of ES cells. 

Gene expression in Rexl-depleted ES cells 

To analyse gene expression in the 7tex7-depleted, AP+ 
colonies, we assayed mRNA for expression of known pluri- 
potency markers or (Polycomb-regulated) differentiation 
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Figure 1. Attenuation of Rexl levels causes partial loss of self-renewal in mouse ES cells. (A) Western blot to detect HA-tagged REX1 in 293T cells 
co-transfected with plasmids that express control shRNAs (con) or two different shRNAs directed against Rexl (shIRexl and sh2Rexl). The amount 
of transfected plasmids is indicated above the panel. A band specific for HA-REX1 is indicated with an arrow, as opposed to several non-specific 
bands. The activity of a co-transfected eGFP expression plasmid served as an indicator of transfection efficiency. The percentage of GFP-positive 
cells is indicated below each lane. (B) Clonogenic colony formation assay. Equal numbers of E14T mES cells transfected with different 
shRNA-expressing plasmids tested in (A) were selected for 7 days for plasmid uptake with hygromycin, and resulting colonies were stained for 
AP. (C) Quantification of colonies obtained in (B). AP+ colonies in a predetermined number of fields were counted; the number of colonies formed 
by cells expressing control shRNA was set at a 100% and the number of colonies in the presence of shRNAs directed against Rexl (shIRexl and 
sh2Rexl) was calculated accordingly. Error bars represent the SD over a minimum of three assays. ** J P<0.01 (ANOVA-Tukey). (D) Changes in 
expression levels of Rexl/Zfp42 after hygromycin-selection for 7 days of cells transfected with plasmids that express shRNAs as described in 
(B). Expression levels were normalized to Gapdh as a reference gene and are depicted as fraction relative to the control ES cells. Error bars represent 
the SD over a minimum of three assays. ***p< 0.001 (ANOVA-Tukey). The difference between effects of shRNAl and shRNA2 is also statistically 
significant, *P<0.05 (ANOVA-Tukey). (E) Western blots to detect REX1 in cell lysates of E14T ES cells (E14T), RA-treated ES cells (RA), or cells 
transfected with shRNAs described in (A). The arrowhead indicates the REX1 protein, which is absent in RA-treated ES cells. Migration of 
molecular weight standards is indicated to the right. Equivalent loading in each lane is demonstrated by stripping and reprobing the membrane 
with aTUBULIN antibodies. 
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Figure 2. Altered expression of ERV upon Rexl depletion in ES cells. Rexl RNA was depleted using transfection of shRNA-expressing plasmids as 
indicated in the legend to Figure IB. After selection for 7 days, changes in gene-expression levels were measured by real-time qPCR: 
(A) pluripotency-associated genes, (B) differentiation markers, (C) several families of ERVs. Expression levels were normalized to Gapdh as a 
reference gene and are depicted as fold expression (2" AACt ) relative to wild-type ES cells (that express control shRNA) (A and C) or fold expression 
(2" AACt ) relative to embryoid bodies (Figure B). Error bars represent the SD over a minimum of three assays. *P<0.05; **P<0.01; ***p< 0.001 
(ANOVA-Tukey). (A) Changes in expression levels of Rexl, Pou5fl/Oct4, Nanog, Sox2 and Tsix. (B) Changes in expression levels of Gata6, Fgf5 and 
T/Brachyury. (C) Changes in expression levels of musD, muERV-L and IAP. 



markers by RT-qPCR. Confirming a previous report (26), 
Tsix RNA was down-regulated to a level proportional to 
the reduction in Rexl levels (Figure 2A). We measured ex- 
pression of several pluripotency markers and detected 
similar or identical levels, independent of the level of Rexl 
for Oct4/Pou5fl and Sox2 mRNAs (Figure 2A). A slight 
increase was observed in the levels of Nanog mRNA upon 
expression of Shi Rexl. This difference was not detected 
upon expression of Sh2Rexl, pointing towards a weak 
dose-dependent effect of Rexl on Nanog expression. No 
up-regulation was observed for several differentiation 
markers previously shown to bind REX1 in ChIP assays: 
Fgf5, Gata6, T (Brachyury) (Figure 2B). We detected small 
differences in expression levels of T upon Rexl depletion 
(Figure 2B), which is consistent with published data 



showing reciprocal expression of Rexl and T (23). We 
conclude that Rexl -depleted ES cells do not exhibit major 
differences in expression of pluripotency markers and dif- 
ferentiation markers. 

RE and ERV are prominently expressed during 
pre-implantation development in the mouse (37,38). As 
Rexl is expressed throughout pre-implantation develop- 
ment (M. Climent et al., submitted for publication) we 
envisaged a potential function for Rexl in regulation of 
ERV. Therefore, we tested de-regulation of ERV in Rexl- 
depleted cells as compared to ES cells with wild-type Rexl 
levels. We initially focused on several autonomous LTR 
retrotransposons, which are highly expressed in either the 
two-cell (muERV-L) or blastocyst stage (IAP, musD) 
mouse pre-implantation embryo. 
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Using primers that amplify regions conserved among 
several copies of REs, we detected a clear difference for 
several families in expression level in Rexl depleted cells as 
opposed to control cells (Figure 2C). Expression of both 
muERV-L and musD was up-regulated ~2- to 3-fold, as 
opposed to IAP, whose overall expression level was not 
affected. We conclude that expression of ERVs is altered 
upon attenuation of Rexl levels in mouse ES cells. 

REX1 binds ERVs of different families to a varying 
degree in mES cells 

We subsequently carried out ChIP assays in combination 
with locus- specific primers to interrogate REX1 associ- 
ation to a subset of ERV in ES cells. Enrichment was 
measured by quantitative real-time PCR and association 
was calculated as percentage association relative to 
control chromatin (Figure 3A). We started out testing 
binding to muERV-L elements, previously demonstrated 
to be highly expressed in the blastocyst (38). After sub- 
traction of background values and normalization against 
an intergenic fragment on chromosome 6 (60), we 
calculated a reproducible 8- to 15-fold enrichment of 
this locus in chromatin immunoprecipitated using 
ocREXl in ES cells (Figure 3 A). To a lesser extent, repro- 
ducible enrichment was also detected for MusD and IAP 
elements (Figure 3A). Similar results were obtained for 
muERV-L and musD in experiments using the Gapdh 
promoter or an intron in the fi-actin gene as controls 
(Supplementary Figure S3). Enrichment for muERV-L 
was also observed when association was compared with 
MLV36 or y-satellite multi-copy sequences, either in 
qPCR (Figure 3B) or semi-quantitative PCR 
(Figure 3C). Weaker, reproducible enrichment was also 
detected for an IAP-related element IAP-RLTR10 
(Supplementary Figure S3). 

Our assay is specific for REX1, as binding was abol- 
ished in ES cells with attenuated levels of Rexl (data not 
shown). Also, we compared binding between undifferenti- 
ated ES cells and their RA-treated counterparts, the latter 
having lost REX1 expression [Figure IE; (27)]. Results 
from ChIP assays (Figure 3D) demonstrate Rexl asso- 
ciation to muERV-L and musD elements in ES cells 
(Figure 3A), as opposed to RA-treated cells (Figure 3D). 

Efficient binding to muERV-L was confirmed by amp- 
lification using different primer sets (Supplementary 
Materials and Methods section), of sequences in the 
5 f -Gag gene or in the 3 f -dUTPase gene (Figure 3E). We 
subsequently tested a potential association of REX1 to 
ORR1 and MT elements of the MaLR family 
(Supplementary Table SIV). While efficient amplification 
was observed using control templates for each of the loci 
analysed (Figure 3E, lanes NoIP), hardly any amplifica- 
tion was observed for the majority of markers after 
immunoprecipitation with pre-immune serum (lane 
Prel). By contrast, we observed reproducible association 
of REX1 to ORR1 and MT elements, as opposed to 
several intergenic regions used as controls (data not 
shown), an intergenic fragment in chromosome 6 
(Figure 3E, bottom panel) and the Gapdh promoter, 
which were not amplified above background (Figure 3E). 



In conclusion, REX1 associates with a subset of 
genomic repeats derived from LTR retrotransposons, 
including muERV-L and IAP elements in mouse ES cells. 

YY1 and YY2 associate with RE 

Although Rexl is expressed in both ES and TS cells in 
culture (27), levels of REX1 association to the same ERVs 
were lower in TS cells as opposed to ES cells (data not 
shown), in line with reduced expression of Rexl in TS cells 
(27). As REX1 shares high homology in the DNA-binding 
zinc fingers with YY1 and YY2 (24), we also tested their 
binding to ERVs. Hardly any enrichment was observed 
for several control markers, i.e. a genomic region in 
chromosome 8 (Genomic A/B), the Gapdh promoter or 
different multicopy controls (MLV36, y-satellites) in 
both cases (Figure 4). After normalization against a 
control promoter, we calculated reproducible enrichment 
of several ERV elements. In chromatin immunopre- 
cipitated using ocYYl (Figure 4A), we demonstrate a 3-, 
5- to 4- and 5-fold enrichment of both IAP and muERV-L 
elements, and a weaker 2-fold enrichment of musD 
elements. We interrogated YY2 association in TS cells, 
considering its high expression in this cell type (data not 
shown). In chromatin immunoprecipitated using ocYY2 
(Figure 4B), we detected a 6- to 8-fold enrichment of 
both IAP and muERV-L elements, and a weaker but sig- 
nificant enrichment of musD elements. In contrast to 
REX1 in ES cells, association of both YY1 and YY2 to 
IAP was slightly stronger when compared with muERV-L, 
indicating that the binding specificities are not the same as 
the one observed for REXl. 

REX1 associates with LSD1 

Transcriptional repression by YY1 is mediated by inter- 
action with HDACs (61). To provide more mechanistic 
insight into the Rexl -dependent repression of retroviral 
elements, we wished to address a potential interplay 
between HDACs and Rexl. We therefore tested the 
effect of the HDAC inhibitor TSA on regulation of 
muERV-L expression. As reported before (46), TSA treat- 
ment caused a significant increase in muERV-L mRNA 
levels in wild- type and Rexl -deficient ES cells (Figure 5A). 
Surprisingly however, the increase caused by TSA was 
attenuated in 7tex7 -deficient ES cells, suggesting HDACs 
function in concert with Rexl (Figure 5A). This effect of 
Rexl was specific, as Rexl deficiency did not influence 
TSA-mediated induction of Dubl (data not shown). The 
fold induction by TSA in Rexl +I+ versus Rexl~^~ cells 
was 1.90 and 1.17 for muERV-L and Dubl, respectively. 
In agreement with recent reports (46), we did not observe 
differences in levels of muERV-L expression in the pres- 
ence of the DNA demethylating agent 5-azacytidine 
(5Aza) (Figure 5A). The resistance to 5Aza was not 
affected by Rexl deficiency, in line with DNA 
methylation-independent pathways being responsible for 
repression. 

Methylation-independent silencing of muERV-L retro- 
viral elements in ES cells is dependent on the lysine 
demethylase LSD1 (46). Other LSD1 target genes are 
also de-regulated in 7tex7-deficient ES cells (22,46). 
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Figure 3. REX1 binds RE/ERV elements in ES cells. (A) REX1 association to RE in ES cells. Binding was assessed in ChIP assays using aREXl 
serum followed by quantification of precipitated DNA using real-time qPCR amplification. The figure shows analysis of a non-binding reference site 
on chromosome 6 as a control (intergenic), and the ERV elements indicated. Association was assessed to the Gag gene in muERV-L (Supplementary 
Materials and Methods section). Association is represented as percentage bound (relative to purified chromatin extract from the same lysate). Error 
bars indicate SEM; *P<0.05. (B) REX1 binding does not extend to other repeated elements. ChIP assays as in (A) except that binding to several 
sequences present in the genome as multiple copies (MLV36, y-satellite) is represented as fold binding relative to a non-binding reference gene Gapdh. 
Error bars indicate SEM. (C) REX1 binding does not extend to other repeated elements. Binding of REX1 in ES cells to several sequences present in 
the genome as multiple copies (MLV36, y-satellite) was assessed by PCR analysis after ChIP by Pre-immune serum (Prel), or aREXl serum 
(aREXl). The figure shows repeat-specific PCR amplification, using primers specific for the elements indicated on the right. PCRs without input 
DNA served as a negative control (MQ). PCRs using a fraction of purified chromatin extract from the same lysate are shown as positive controls on 
the left (No IP). (D) Loss of REX1 binding to selected ERV elements in RA-treated ES cells. Levels of REX1 binding were determined by ChIP as in 
A using either pluripotent or RA-treated ES cells and site-specific qPCR. The relative amount of immunoprecipitated DNA is shown as the ratio of 
aREXl immunoprecipitations in RA-treated versus non-treated cells, normalized to non-bound control genes (Intergenic or Gapdh). Error bars 
indicate SD. (E) REX1 association to muERV-L and MaLR. Binding of REX1 in ES cells was assessed by PCR analysis after ChIP by Pre-immune 
serum (Prel), or aREXl serum (aREX). The figure shows gene-specific PCR amplification, using primers specific for the ERV indicated on the right. 
PCRs without input DNA served as a negative control (MQ). PCRs using a fraction of purified chromatin extract from the same lysate are shown as 
positive controls on the left (NolP). 
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Figure 4. YY1 and YY2 association to ERV elements. Association of YY1 in ES cells (A) and YY2 in TS cells (B) to potential target genes was 
assessed by /ocws-specific qPCR analysis after ChIP using aYYl or aYY2, respectively. The Gapdh promoter, non-binding control elements in 
chromosome 8 and several sequences present in the genome as multiple copies (MLV36, y-satellite) are included as negative controls. Association was 
assayed for musD, IAP and muERV-L. The amount of immunoprecipitated DNA as a percentage of input DNA was recalculated as fold association 
normalized to the Gapdh promoter. Error bars represent SEM over a minimum of three assays. 



As these data suggest co-operation between REX1 and 
LSD1, we probed for direct interactions. To do so, we 
introduced plasmids that drive expression of either HA- 
tagged REX1, together with plasmids that express 
cDNAs encoding LSD1 protein (48) in tissue culture cells. 
Interacting proteins were identified by western blot after 
immunoprecipitation. Using ocREXl serum, we detected 
LSD1 co-immunoprecipitated with HA-REX1 (Figure 5B, 
lanes marked LSD1). Detection of the interaction was de- 
pendent on the use of ocREXl serum as opposed to 
Pre-immune serum (Supplementary Figure S5B), and on 
the presence of REX1 proteins (Figure 5B). Furthermore, 
in the reverse experiment, we easily detected HA- REX 
associated with immunoprecipitated Flag-tagged LSD1 



(Supplementary Figure S5C). These results indicate that 
REX1 is capable of association with LSD1 protein. We 
propose that REX1 may function in regulation of 
muERV-L expression in concert with LSD1 and associated 
HDACs. 

Rexl represses expression of ERV elements in gain- and 
loss-of-function mouse embryos 

As we had observed an increase in expression of LTR 
retrotransposons in 7tex7 -depleted ES cells when 
compared with wild-type, we measured directly whether 
de-repression of retroelements was associated with mobil- 
ization (Supplementary Figure S4). As Rexl participates 
in Polycomb regulation and on the other hand Rnf2 
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Figure 5. Association of REX1 and LSD1 proteins. (A) Expression 
levels of muERV-L. Rexl +I+ or Rexl~^~ ES cells were treated with 
TSA or 5Aza for 24 h. Levels of muERV-L expression were determined 
by qRT-PCR as explained in the legend to Figure 2C. Results are 
depicted as fold expression in each cell line after normalization to 
Gapdh. The fold change was then plotted relative to untreated cells, 
with error bars representing SD. The fold muERV-L induction by TSA 
in Rexl + ^ + versus Rexl~ ! ~CQ\h was 1.90. (B) Human kidney 293T cells 
were co-transfected with a plasmid expressing HA-tagged-intact 
REX1 proteins as indicated, in combination with a plasmid expressing 
LSDl/Kdmla. Proteins were immunoprecipitated with aREXl 
antibodies and analysed by western blot (top panels) together with 
proteins in total extracts (bottom panels) with monoclonal aHA 
antibody, aREXl antibodies or aLSDl antibodies as indicated on 
the right. Bands corresponding to HA-REX1 or LSD1 are marked 
on the right, protein molecular weight markers are indicated to the 
left in the bottom panels. 



deletion was reported to severely affect retro-transposition 
of MLV36, a mouse leukaemia virus (ERV Class I) (62), 
we also tested mobilization of this retrovirus in Rexl- 
depleted ES cells. 

We did not observe increases in the copy number 
of IAP, musD, Etn, muERV-L elements or of the 
MLV provirus in ES cells with attenuated levels of 
Rexl as opposed to control cells (Supplementary Figure 
S4). We infer that although musD and muERV-L 
elements were overexpressed when Rexl was depleted 
(Figure 2), they could not integrate into the genome at 
increased levels. 

MuERV-L is expressed during pre-implantation devel- 
opment and displays a transient increased expression 
during the two-cell stage (38). To prove a potential role 
for Rexl in control of muERV-L expression, we generated 
Rexl loss-of-function embryos by micro-injection of 
zygotes with siRNAs directed against Rexl. Embryos 
were further cultured in vitro, and gene expression was 
assayed by RT-PCR at different stages using expression 
of H2afz as a reference gene, comparing embryos injected 
with control siRNAs or anti-7tex7 siRNAs (see 'Materials 
and Methods' section). As opposed to H2afz, expression 
of Rexl was clearly affected by the siRNAs (Figure 6A, 
lanes 3 and 5). Rexl mRNA levels were moderately 
attenuated in both two-cell embryos and morulas. 
Under these conditions, muERV-L levels were not 
altered in two-cell embryos (Figure 6 A, lane 3). By 
contrast, we observed a significant difference in morula, 
as Rexl attenuation resulted in ectopic expression of 
muERV-L. 

To test the effects of Rexl gain-of-function, we also 
injected zygotes with expression vectors that overexpress 
Rexl or eGFP as a control (See 'Materials and Methods' 
section). After continued development in vitro, we selected 
embryos with either high or low expression of GFP (from 
an attached IRES-eGFP in the case of Rexl) and assayed 
differences in gene expression between Rexl and eGFP 
overexpressing embryos by RT-PCR. No differences 
were observed between wild-type, eGFP overexpressing 
embryos and embryos that express absent or low levels 
of the injected 7tex7-IRES-eGFP plasmid (data not 
shown). Therefore, we compared embryos that express 
either neglectable or high levels of eGFP, after injection 
with 7tex/-IRES-eGFP plasmid. Results are presented in 
Figure 6B. We easily observed expression in blastocysts of 
H2afz (control), IAP and muERV-L as described (38). 
Overexpression of Rexl in these embryos (Figure 6B, 
lane 3) suppressed expression of both IAP and 
muERV-L elements. This reduction was not caused by a 
general effect on transcription, as expression of Oct4\ 
Pou5fl was not affected by Rexl overexpression. 

We conclude that attenuation of Rexl levels interferes 
with repression of ERVs, while overexpression attenuates 
normal expression levels. These results demonstrate that 
Rexl controls expression of ERV, particularly of the 
muERV-L family, during pre-implantation mouse devel- 
opment. These data indicate that Rexl has a crucial role in 
controlling ERVs both during early embryonic develop- 
ment and in ES cells. 
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Figure 6. Rexl gain-and-loss of function in mouse embryos. (A) Effect of attenuation of Rexl levels on muERV-L expression. One-cell embryos 
were microinjected with control siRNAs (C) or 7tex7 -specific siRNAs (R), cultured to either the two-cell or the morula stage and analysed by 
semi-quantitative RT-PCR. cDNA from equivalent amounts of embryos was amplified with primers specific for the genes indicated; H2afz was used 
as a control for input. The top panel shows an ethidium bromide stained gel of amplification products. No I refers to no cDNA input. RNA 
extracted from E14T ES cells was processed alongside as a positive control. Bottom panel: quantification of Rexl and muERV-L levels (grey and 
black bars for C and R, respectively). Data are represented as the average fold induction in Rexl depleted morulas (black bars) using H2afz as a 
reference gene (two independent experiments). (B) Top panel: muERV-L and IAP mRNA expression in wild-type or Rexl overexpressing mouse 
blastocysts. One-cell embryos were microinjected with plasmids that direct expression of 7tex7-IRES-eGFP. Embryos were separated in groups 
according to eGFP expression levels at the eight-cell stage, and 4 days after injection RNA expression was analysed by semiquantitative RT-PCR as 
described in Figure 6(A). The following templates were used for amplification: E14T ES cells (lane 1), embryos with low or neglectable Rexl 
overexpression levels (lane 2), embryos with high Rexl overexpression (lane 3) or MQ water as a negative control (lane 4). Bottom panel: quanti- 
fication of Rexl, Pou5fl/Oct4 and muERV-L levels, represented as fold repression in Rexl overexpressing embryos using H2afz as a reference gene. 



DISCUSSION 

Rexl regulates expression of ERVs 

We show that attenuation of Rexl levels de-regulates ex- 
pression of ERV elements both in mouse ES cells and 
embryos. As REX1 also specifically binds these elements 
in mouse ES cells, we provide evidence for regulation by 
Rexl of ERVs, especially muERV-L elements. To a lesser 
extent, binding was also observed to musD, IAP elements. 
It should be taken into account that the primers used were 
designed against conserved regions and therefore amplify 
several copies of each family. We suggest that this 
property may underestimate changes in both gene expres- 
sion and association of REX1, as non-regulated copies 
may also be assayed. In mouse ES cells and in the presence 
of attenuated but detectable levels of Rexl (Figure 1), 



regulation did not involve retro-transposition itself 
(Supplementary Figure S4). We cannot rule out the pos- 
sibility that during pre-implantation development, or in 
ES cell truly deficient in Rexl, retro-transposition is 
affected by Rexl. 

We demonstrate an important role for Rexl in control 
of ERV, as Rexl depletion in pre-implantation embryos 
provoked up-regulation of muERV-L, while the gain-of- 
function experiments caused down-regulation of several 
ERVs, including muERV-L (Figure 6). The de-regulation 
of muERV-L (and IAP) ERVs represents the earliest 
phenotypic abnormality described in Rexl loss of 
function embryos. 

During pre-implantation development, REs and par- 
ticularly ERVs are prominently expressed in a precise 
temporal pattern (37-40), while general demethylation 
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(41) removes the major repressing mechanisms active in 
adult tissues (63). Control of ERV expression during 
pre-implantation development involves protection of 
methylation patterns imposed by Dppa3j Stella (44) and 
the expression of ERV-specific siRNAs (43). 

Recently, KDM1 A/LSD 1 has been discovered as a 
methylation-independent epigenetic regulator of ERV 
and ERV LTR-linked genes (46). Reminiscent of Rexl, 
Ls<i7-deficient ES cells could be derived and deficiency 
caused de-regulation of pre-implantation-specific gene ex- 
pression, in turn interfering with lineage determination 
upon differentiation. In addition, genome-wide analysis 
of binding sites in ES cells has revealed preferential 
binding of REX1 to RE-linked genes expressed during 
pre-implantation development (Guallar,D., Garcia- 
TunonJ., Climent,M., Muniesa,P. and SchoorlemmerJ. 
manuscript in preparation). LSD1 is associated with 
HDACs (46,64) and HDAC-mediated repression of 
muERV-L retroviruses in ES cells is affected by Rexl de- 
pletion (Figure 5A). Although we show that REX1 and 
LSD1 can interact (Figure 5B), the interrelationship of 
REX1 with either of these mechanisms in the control of 
ERV expression has not been molecularly defined yet. We 
speculate that similar to or even in concert with LSD1 
and/or associated HDACs, REX1 may play a general 
role in control of transcription from ERV and ERV- 
derived elements during mouse pre-implantation 
development. 

Rexl and self-renewal/pluripotency in ES cells 

Although it has been reported that Rexl null mice appear 
normal, viable and fertile and that normal ES cells can 
readily be derived either from Rexl null mice or by hom- 
ologous recombination in vitro (22,23), an important role 
for Rexl in ES cells is suggested by the high coincidence of 
Rexl expression and enhanced pluripotency (18,20). We 
observe that the attenuation of Rexl levels negatively 
affects self-renewal in mouse ES cells (Figure 1), in 
contrast to results reported earlier (22,23). The discrep- 
ancy may be due to functional redundancy with Yyl 
and Yy2 proteins. Alternatively, the diminished self- 
renewal ability caused by depletion of Rexl levels we dem- 
onstrate may only be revealed in the clonogenic assays we 
use. Surprisingly, we found no evidence that expression of 
muERV-L is elevated in Rexl -deficient ES cells (data not 
shown), and TSA treatment relieved repression of 
muERV-L in Rexl -deficient cells (Figure 5A). Therefore, 
in addition to repression through a 7tex7/HDAC- 
dependent mechanism, a parallel pathway is operative to 
repress muERV-L in mouse ES cells. As repression of 
muERV-L is not notably influenced by demethylating 
agents, even in the absence of Rexl, we suggest that in 
ES cells neither mechanism is dependent on DNA methy- 
lation. We propose that a potential player in the second 
mechanisms might be RYBP (65). 

As we fail to detect major differences in ES cell-specific 
gene expression upon Rexl depletion, we suggest that, 
rather than directing pluripotency- and self-renewal- 
related gene expression, Rexl may contribute to 
pluripotency and pre-implantation development via 



mechanisms that affect or read chromatin structure 
around ERVs. Such a role is also exemplified by the regu- 
lation of imprinted genes by Rexl (28) at stages preceding 
the generation of pluripotent cells and by the epigenetic 
role of Rexl in de-regulation of differentiation (22), 
possibly through association of REX1 to Polycomb 
group proteins (27). 

ERVs and retroviral LTRs have been co-opted by 
cellular genes as promoter elements. In this way, ERVs 
contribute to -25% of POU5F1/OCT4 and NANOG- 
binding sites, both in human and mouse ES cells (66). 
ERVs, and especially ERV-K (IAP, musD, ETn), are 
over-represented among binding sites for NANOG, 
SOX2 and OCT4 (67). We speculate that Rexl interaction 
with the principal pluripotency network (10) is centered at 
RE/ERVs and/or LTR-linked genes. Indirectly, Rexl may 
influence posterior differentiation, similar to endo- 
dermal differentiation resulting from KDM1 A/LSD 1 
deficiency (46). 

A function for YYl-family members 
in control of REs/ERVs? 

We show that REX1 is significantly enriched at genomic 
loci encoding members of the ERV2 and ERV3 
superfamilies of LTR-containing ERVs) in mouse ES 
cells (Figure 3 and Supplementary Table SIV). Specifi- 
cally, we observed binding of REX1 to muERV-L, 
ORR1 and MT elements, and to a lesser extent to IAP 
and musD elements. We do not know at present to which 
extent REX1 associates with a wider subset of ERV, but 
conclude that REX1 associates with at least a distinct 
subset of ERVs. Searches for REX 1 -binding sites within 
IAP/MusD/muERV-L sequences using position weight 
matrices revealed the presence of several potential 
binding sites, suggesting that regulation likely results 
from direct association of REX1 to DNA (B. Moreiras, 
data not shown). More extended analysis will be required 
to determine whether these associations are (super)family- 
dependent or are determined by the sequence and/or epi- 
genetic environment of individual elements. Alternatively, 
the association of REX1 may be tissue-specific, as relative 
affinity for IAP and muERV-L elements is changed in 
placental tissue (M. Climent, R. Perez-Palacios, 
S. Climent, P. Muniesa and J. Schoorlemmer, unpublished 
data). 

REX1 has been described to activate Tsix expression 
through binding a ds-acting element discovered as 
DxPas34 (68,69). Not surprisingly, DxPas34 shares 
homology to the ERV3/ERV-L family of ERVs (69,70). 
Homology to this element is conserved in and restricted to 
human, mus musculus, rodents and eutherian mammals in 
general. These homologies support an origin of DxPas34 
as an ERV-originated element and are compatible with a 
generic function for Rexl in control of transcription from 
c/s-elements that originate from ERVs. It has been sug- 
gested that ERV-L inserted more than 70 million years 
ago in the Tsix locus and was adopted as a regulatory 
module for Xi by early eutherians (69). Interestingly, 
Rexl was duplicated from Yyl by retro-transposition. 
This event also happened exclusively in eutherian 
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mammals about 60-100 million years ago (24). We 
propose that similar to co-option of ERV as a regulatory 
element for Xi, Rexl was co-recruited as a regulator of 
such elements. 

All three Yyl family members are present during 
pre-implantation development [(13,71), R. Perez and J. 
Schoorlemmer, unpublished data] and in the germ line 
(13,72). REX1 shares extensive homology in the 
DNA-binding zinc fingers with YY1 and YY2, which 
prompted us to test association of YY1 and YY2 to 
ERVs bound by REX1 (Figure 4). We demonstrate ap- 
preciable binding of both proteins to muERV-L and IAP 
elements, hinting at a shared function of family members. 
As specificity was slightly altered compared to REX1, we 
hypothesize that each protein may regulate a specific 
subset of REs. De novo binding sites generated by trans- 
posons have contributed significantly to species-specific 
gene regulatory networks during pre-implantation devel- 
opment (67). As ERVs are less active in human as opposed 
to mouse (31) and certainly less prominent in the region 
corresponding to DXPas34 (69,70), it will be of interest to 
survey the RE specificity of YY1 -family members in 
human. 
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